On the description of exclusive processes beyond the leading twist approximation 
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We describe hard exclusive processes beyond the leading twist approximation in a framework based 
on the Taylor expansion of the amplitude around the dominant light-cone directions. This naturally 
introduces an appropriate set of non-perturbative correlators whose number is minimalized after 
taking into account QCD equations of motion and the invariance under rotation on the light-cone. 
We examplify this method at the twist 3 level and show that the coordinate and momentum space 
descriptions are fully equivalent. 
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Exclusive reactions have been the scene of much re- 
cent progress in the domain of perturbative QCD physics. 
CoUinear factorization of the amplitude for vector meson 
production has been demonstrated at the leading twist 
level [0. In this case the meson field appears through 
twist 2 distribution amplitudes (DAs) which are defined 
as the Fourier transforms of the matrix element of nonlo- 
cal two particle operators 0{x), such as tlj{x)Tip{Q) for 
fermionic partons, F being an appropriate Dirac struc- 
ture. Generalizing this factorization procedure to the 
twist 3 case is mandatory in a number of cases (provided 
the value of the hard scale is sufficiently high to justify 
phcnomcnologically the twist expansion), and in partic- 
ular for the understanding of the production of transver- 
sally polarized vector mesons ||^ , since their twist 2 DAs 
are often decoupled because of their chiral-odd nature [|j . 

Following the strategy outlined in Q|, we here gener- 
alize the approaches ^ in order to get a thorough 
description of the exclusive processes. In particular, wc 
discuss the role of the light-cone direction independency 
condition and fully explore its consequences in the anal- 
ysis of genuine higher twist effects, in particular beyond 
the Wandzura Wilczek approximation. We show that 
the present framework is equivalent to the coordinate de- 
scription developed in [|[ Q . 

In this paper we summarize only the main ideas of our 
developments, which are to a big extent valid for arbi- 
trary exclusive processes. However, for illustration of the 
main steps of the derivation, in particular of the imple- 
mentation of the invariance under rotation on the light- 
cone of the physical amplitude, we will rely on the specific 
meson impact factor case at twist 3, when needed. This 
quantity is illustrated in Fig. |l|, it enters the descrip- 
tion of transversly polarized vector mesons production 
at small x, i.e. at large energy, in 7*7* and j*p col- 
lisions. For the aims of present paper it is enough to 
understand that it represents a part of the amplitudes of 
above mentioned processes which describes the transition 



of a virtual photon to a vector meson mediated by two 
virtual gluons being in the colorless state. Full details 
of the derivations and explicit definition of impact factor 
arc given in |^ . 

Thus, we start with the amplitude presented in Figj^, 
written as (with {d^£}2 = d'^h ^^^2) 



A= d'^i tr 



id' 



.tr 



(1) 



where H and iJ^ arc the hard parts with two parton 
legs and two parton and one gluon legs, respectively. It 
is important to note that it is assumed that virtualities 
of i— channel gluons in Fig. |l| are large, being of order 
of photon virtuality Q^. Thus, all nonperturbative, soft 
momenta part of the subprocess described by the impact 
factor is included in the blobs $ and In (|^), the soft 




FIG. 1: 2- and 3-body correlators attached to a hard scatter- 
ing amplitude in the example of the 7* ^ p impact factor. 

parts are given by the Fourier-transformed two or three 
body correlators. Inclusion of 4-body correlators would 
lead to twist > 4 contributions, which are suppressed by 
additional powers of 1/Q. Besides, here wc consider the 
case of vanishing transverse momentum transfer in the 
7* — > p transition and we neglect corrections of the order 
rUp/Q. 

Wc have to parameterize momenta and polarization 
vector of the transversly polarized meson (e). In our 
subprocess 7*((j) p{Pp), = 



Q2 



Pp^ 



(2) 
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Here p (n) be a "plus" ("minus") light-cone vector, re- 
spectively (normalized as p ■ n = 1). Such light-cone 
vectors can be easily constructed from the 4- vectors de- 
scribing the kinematics of the process in question. 

The polarization vector of transversly polarized meson 
obeys 



e ■ p ~ Q 



e • n 







(3) 



This vector coincides with the physical polarization vec- 
tor in the infinite p— meson momentum frame. 

The amplitude (|l|) is not factorized yet because the 
hard and soft parts are related by the 4-dimensional in- 
tegration in the momentum space and by the summa- 
tion over the Dirac and/or Lorentz indices. To factorize 
the amplitude, we need to select the dominant direction 
around which we intend to decompose our relevant mo- 
menta and the hard part. For this purpose we choose 
another pair of light-cone vectors: p and n (normalized 
again as, p ■ n = \). It is very important to note that 
here we take the freedom to choose n arbitrary, having 
no direct relationship with the parametrization of exter- 
nal momenta, n ^ fi. We carry out an expansion of £ in 
the basis defined by p and n: 



VtPt, + i^t -p) n^,+ii 



(4) 



with Hi ~ ii ■ n , and rewrite the integration measure in 



(1^ as dH,, 



(rlidyi5{yi 



We also decompose 



the hard part around the dominant "plus" direction: 



H{i) = H{yp) 



dH{e) 



(5) 



with {£ ~ yp)a ~ £a up to twist 3 accuracy, where a 
denotes transverse components with respect to vectors p 
and n. 

The dependence of the hard part looks first as an 
excursion out of the coUinear framework. The factorized 
expression is obtained as the result of an integration by 
parts which replaces by acting on the soft cor- 
relator. This leads to new operators which contain 
transverse derivatives, such as ip d'^ip, and thus to the ne- 
cessity of considering additional correlators ^'^{l)- This 
procedure leads to the factorization of the amplitude in 
momentum space. Factorization in the Dirac space can 
be achieved by the Fierz decomposition. For example, 
in the case of two fermions, one should project out the 
Dirac matrix ipa{0)ip[3{z) which appears in the soft part 
of the amplitude on the relevant F matrices. Thus, the 
amplitude finally takes the simple factorized form: 

1 

A = J dytr[Hiy)T] (y) 



1 

+ J d2/irfy2tr[ff^(yi,2/2)r] ^l{yi,y2)- (6) 



The first term of Eq. (|) requires the introduction of the 
soft correlators with quark-antiquark nonlocal operators. 
Generically jiot , 

{pipMz)j^^pmo) 



= mp fp [Viiy) (e* • n)pf, + ip^ijj) e* ] , (7) 



(p(p)|V^(z)757pV^(0)|0) 



mpfpiipA{y)ef,xpsel'^Pl3ns, (8) 



:Fi 



where — denotes the Fourier transformation 
/p dyexp[iyp-z] . The light-like separation z ^ \n 
involves an arbitrary light-like vector n with n ■ p = 1 
which already appeared in the Sudakov decomposition 
(^). The above correlators are defined in the axial 
light-like gauge n ■ A = , which allows to get rid of 
Wilson lines. Thus, n plays a triple role: it first fixes the 
axial gauge, second it defines the transverse momentum 
of partons forming the p in the coUinear expansion, 
and last it determines the notion e/^ of transverse 
polarization of the p in the above equations. Note that 
should not be confused with the physical polarization 
vector of transverse meson, defined by the conditions 

The second term of Eq. (|^) requires the introduction 
of matrix elements of quark-transverse gluon operators 

(9) 



{p{p)\4'iz,)-,^gAi{z2)mm 

•mpf3pB{yi,y2)p,j.e^ 

{pipMizih5l^.9Al{z2)m\0) = 



a ' 



(10) 



■mp hpiD{yi,y2)pp,ea\fise*x ppns, 

1 1 

where / dyi J dy2 exp [iyip- zi + i{y2 -yi)p- Z2] is de- 


noted by = . Note that the positivity of the gluon light- 
cone momentum fraction imposes that quark-giuon pa- 
rameterizing functions have the form 

B/D{yi,y2) = B/V{yi,y2;y2 - yi) d{yi < y2 < 1) ■ 

Here, the light-cone fractions of the quark, anti-quark 
and gluon are respectively yi , 1 — ^2 and j/2 ^ J/i • 

Finally, the second term of (H) leads to the introduction 
of the correlators 

{p{p)mz)jpi d;t mm = mpfp^f{y)p^e*J ,{U) 

{p{pmzh5lp^ mm 

= mpfpiif>'^{y)p^ea\f3se*yjpf3ns. (12) 

Let us recall the symmetry properties of DAs Q com- 
ing from C invariance (with y = 1 — y): 

vi{y) = fiiy) , f3{y) = 'P3{y), fAiy) = -fAiy) , 
v^f (y) = -vf (y) , VAiy) = y^A{y), 

B{yi,y2;yg) = ~B{y2,yi;yg) , 

V{yi,y2;yg) = V{y2,yi;yg). (13) 
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The correlators introduced above are not independent, 
thanks to the QCD equations of motion for the field 
operators entering them (see, for example, Q]). In the 
simplest case of fermionic fields, they follow from the 
vanishing matrix elements {{iD{0)ip{0))a '4)p{z)) = and 
{•ipa{^)'i{D(z)'tp{z))p) = due to the Dirac equation, 
then projected on different Fierz structure. Denoting 



V,A 
3,p 



f^flfp^ obtain 



'3p 



yi fsivi) + yi VA{yi) + v^iyi) + VAiyi) 
1 

^- j '^y^ [^^ ^(yi' 2/2) + Ct D{yi, y2)] , (14) 



yi miyi) ~ yi 'PAiyi) - (yi) + VaIvi) 
1 

= - I dy2 [-Q^ B{y2, yi) + C3 D{y2, yi)] ■ (15) 



Note that Eq.(|l5|) can be obtained by the replacement 
yi — > yi hi ( [1^ ) and the use of symmetry properties (13). 

We stress again that contrarily to the light-cone vector 
p related to the out-going meson momentum, the second 
light-cone vector ri, required for the parametrization of 
the correlators (0-^2|), is completely arbitrary and con- 
nected in no way to the momenta of particles of the 
process in question. Alternatively, the notions of trans- 
verse momentum, of transverse polarization and of axial 
gauge, lead to a description of the factorization proce- 
dure which is non-covariant. One should thus impose, 
due to Lorentz invariance, that any physical observables 
do not depend on the specific choice of n, e.g. the scat- 
tering amplitudes should be n— independent. It turns 
out that this requirement leads to additional non-trivial 
constraints between the non-perturbative correlators en- 
tering the factorized amplitude. This independency 
property was crucial for the understanding of inclusive 
structure functions properties at the twist three level ||] 
and its relevance for some exclusive processes was pointed 
out in [Q. We show now that this condition expressed 
at the level of the full amplitude of any process can be 
reduced to a set of conditions involving only the soft cor- 
relators. The strategy relies on the power of the Ward 
identities to relate firstly amplitudes with different num- 
ber of legs and secondly higher order coefficients in the 
Taylor expansion (|^) to lower order ones. 

At first sight the n— independency condition of the pro- 
cess amplitude could not lead to some process indepen- 
dent results. Indeed, not only the parameterizations of 
the hadronic matrix elements but also the process depen- 
dent hard scattering coefficient H depends (indirectly) on 
n, since n defines the direction of the coUinear expansion, 
see eq. (||). With this respect we remind that using in- 
tegration by parts in the d^l integral (|^) one replaces 
by acting on the soft correlator. Therefore the n— 



independency of the process amplitude can be reduced 
entirely in considered approximation (at twist-3) to the 
consideration of the n— dependency of the soft correla- 
tors, see eqs. (p|j7 Hll|JT^ ,p| JTo|), which parameterizations 
depend on n both explicitly and implicitly through the 
n— dependence of the polarization vector e^. 

The light-cone vector n satisfying n ■ p = I may be 
parametrised by its transverse components n± defined 
with respect to the light-cone vectors p and h fixed by 
the external kinematics (see eq. (pT)) 



■P 



(16) 



Thus the n— independence of A for an arbitrary fixed 
polarization vector e is expressed by the condition 



(17) 



On the other hand, the scattering amplitude A receives 
contributions from the vector correlators, which result 
into ^'"'2c*or^ from the axial vector correlators, which 
lead to ^f*^*"^ part of A. Due to different parity prop- 
erties of the vector and the axial-vector correlators, the 
condition (fi^) means effectively two separate conditions: 



and 



dn'1 



A" 



A"" 











(18) 



(19) 



The dependence of A on the vector n± is obtained 
through the dependence of A on the full vector n. This 
dependence on n is different in A"'^'''^^ and in _4'"^'^*°'' 
parts. 

The dependence of A"'^^"''' on the vector n enters only 
through the expression £P"'^t involving the contraction 
with the momentum p and in which the indices /? and 7 
are contracted with some other vectors. Thus the condi- 
tion (|l^) is equivalent to 



d 



dill dn°' 



A 



axial 



j^axial ^ g ^20) 



where we took into account the peculiar dependence of 
jaxiai ^ discussed above. This will lead at the level 
of DAs to the equation 



d 



dy2 



dyi J 2/2-2/1 



X {D{yi,y2) + D{y2,yi)) . 



(21) 
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(j/i - 2/2 ) 



FIG. 2: Reduction of 3-body correlators to 2-body correlators 
through Ward identity. 



The polarization vector for transverse p which enters in 
the parametrization of twist 3 correlators is defined as 



= - e • n . 



(22) 



It turns out, as one can check explicitely in this example, 
that the dependence of on the vector n enters 

only through the scalar product e* • n, and Eq. ( |l8| ) can 
be written as 



(7 r) 

-vector g* M 



A'' 



d (e* • n) 



A" 



= 



(23) 



from which follows 
d 







dy2 



X {B{yi,y2) + B{y2,yi)) 



(24) 



We now sketch the derivation of Eq.(gJ), which cor- 
responds to the vector correlator contributions, in the 
peculiar case of the 7* — > p impact factor, which is the 
building block of the description of the j* p — > pp pro- 
cess at large s . The 3-body {qqg) contribution and the 
2-body contribution involving to A can be reduced to 
the convolution of the leading order hard 2 particle con- 
tributions with linear combination of correlators, thanks 
to the use of the Ward identity. In the case of the 3-body 
vector correlator (||), due to ( p^ ) the dependency of er 
on n_L enters linearly and only through the scalar prod- 
uct e* • n . Using the fact that the amplitude A depends 
linearly on the polarisation vector e , one finds that the 



action on the amplitude of the derivative 
can be extracted by the replacement e* 
using the Ward identity, it reads 



d 



d (e* -n) 



-Pa ■ After 



(2/1 



-y2)tr[Hp{yi,y2)pP^]=tr[H{yi)^]-tr[H{y2)^] , 

as illustrated by Fig.^. This results in the last term of 
r.h.s of Eq.(p^). A similar treatment of 2-body correla- 
tors with transverse derivative whose contribution can be 
viewed as 3-body processes with vanishing gluon momen- 
tum leads to the l.h.s of Eq.(^J). The first term with (pi 
of the r.h.s of Eq.(^) originates from the 2-body vector 
correlator and corresponds to the contribution of ^ P- 
The second term with p^ corresponds to the contribution 



of the same correlator for pT supplemented by Eq.(22) 



A similar treatment for axial correlators leads to Eq.(21) 



Since we rely on the n— independency of the amplitude, 
one may wonder about the effect of the gauge choice. 



which is fixed by n, on the hard part. The QCD Ward 
identities require the vanishing of the amplitude in which 
polarization vector of a gluon is replaced by its momen- 
tum provided all other partons are on the mass shell. 
In the framework of the fcr— factorization the i— channel 
gluons are off the mass-shell. Therefore the replacement 
of the s— channel gluon polarization vector by its momen- 
tum leads to the vanishing of scattering amplitude up to 
terms proportional to k'^/s where k± are transverse mo- 
menta of t— channel gluons. From the point of view of 
the i— channel, the gauge invariance of the impact-factor 
means that it should vanish when the transverse momen- 
tum of any channel gluon vanishes. To acchieve this 
property it is necessary to include in a consistent way not 
only DAs with lowest Fock state containing only quarks 
but also those involving quarks and gluon, see 

In practice, we here check this invariance by contract- 
ing the channel emitted gluon vertex in the hard part 
with the momentum, which in coUinear factorization is 
proportional to the p— meson momentum, this leads to 
simplifications in the use of (coUinear) Ward identities. 
This means that the hard part is gauge invariant. 

In order to prove this, one should first project on the 
various color Casimir structure. In the case of the impact 
factor, this means to distinguish Nc and Cp terms. In 
this case, Cp terms arise from 2-partons diagrams and 
from 3-partons diagrams where the emitted gluon is at- 
tached to a quark line, while Nc terms are obtained from 
3-partons diagrams where the emitted gluon is attached 
to a quark line only between 2 channel exchanged glu- 
ons or from diagrams involving at least one triple gluon 
vertex. Using Ward identity, the result looks like the 
product of the 2-parton diagrams contributing to the 
transition —f p^ . One can readily check at Born oder 
that this transition vanishes, in our twist 3 treatment. 
In the case of contribution proportional to Nc, one can 
prove that this gives zero. Note that terms proportional 
to the virtuality k^ remains, but they are of the form 
k^/s , which are higher order with respect to the domi- 
nant s— power contributions to the impact factor (which 
scales like s'^) involved in the fcy— factorization approach. 
This is thus in agreement with Ward identity in pertur- 
bative Regge factorization. 

Although our implementation of factorization and 
n— independence condition was illustrated on the partic- 
ular example of the impact factor at twist 3, we expect 
that this procedure is more general and that the above 
method can be applied for other exclusive processes, the 
key tool being the coUinear Ward identity. This means in 
particular that each building block (soft and hard part, 
for each structure which lead to the introduction of a 
DA) are separately gauge invariant. This fact simplifies 
dramatically the use of the rt— independence principle. 

We thus emphasize that the equations obtained above 
for the special case of impact factor are expected to 
be universal i.e. not to refer to any specific hard pro- 
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cess. The equations of motion (^,|T^) together with Eqs. 
(^,|2j) following from n-indcpendcncc arc four equations 
which constrain seven DAs (/^a, (p^, (^^, (p^, B, D). 
Thus, any hard process involving the exclusive produc- 
tion of a p can be expressed, at twist 3, in terms of the 
three independent DAs which we choose as </ji, B and 



D. The solution of Eqs. (y,|l5|,|24|^ with vanishing 
B and D DAs corresponds to the Wandzura-Wilczck ap- 
proximation and takes the form |12| 



d d 
I —fi{v) ± / 

J V J V 

.0 y 

y 1 



y 

The remaining genuine twist 3 DAs read: 



(25) 
.(26) 



gen / \ 



y 
1 

{CXB-(:iD){v2.u) 



y2-u 
dy2 



y2~u 



yi 



u 

1 

dy2 



{C^B + C^D){u, y2) 
iC^B + CiD){y2,u) 



T qen / \ 

<Pi yy) 



y2-u 
dy2 

y2-u 



du Lp% (u) - Cs / dyi / dy; 



(27) 



Bjyi, y2) 
2/2-2/1 



(28) 

while the corresponding expressions for ip^^"'{y) and 
'PA^^^iy) ^^'■^ obtained by the substitutions: 



"Pa (y) ' — ' <^3 (2/): 



(29) 



^^^"(2/) ^-^"^^^ ^^^^"(2/). (30) 

Let us now discuss an equivalent way for obtaining the 
result (|2^). We start with the following operator equa- 
tion written in a gauge-invariant form, where on the 
right-hand side wc only keep the contributions of vector 



correlators 



d 

dZa 




= f31) 


~1p{z)^fj,[z, ~z]Da Tp{-Z) 
1 


+ ip{z)Da 7^[z, -z]ip{-z) 




f dv V iIj{z)[z, vz\zvG 


va{vz)"f^[vz, -Z\'4){-Z) , 




1 





where [] denotes the Wilson line and Da = da —igAa{—z) 

and Da=da +igAa{z). Note that in (|3l|), the derivatives 
act on the arguments of the fcrmion fields. In princi- 
ple, the path in the Wilson line can be an arbitrary one, 
provided it lies entirely on the light-cone. However, for 
simplicity, we choose the path to be a straight line. We 
now compute the matrix elements of the various part of 
Eq.dl). 

First, using the parametrization (^, the matrix ele- 
ment of the quark-gluon part on r.h.s of Eq.(^) reads 

{p{p)\4iz)j^{[z, ~z]Ao,i-z) + A^iz)[z, -z]}^(-z)|0) 



dvv {p{p)\iljiz)[z, vz\zyGua{vz)^^_,[vz, -z]-0(-z)|O) 



2 mpfl„p, e*J dy e'fe-S)f-/ dx, / dx^ 



B{xi,X2) 
X2 - Xi 



(32) 



where in the last stage integration by parts has been used. 

Further, we turn on the quark-antiquark operator in 
( pl| ) and its vacuum-to-meson correlator. Using ([2^), 
the correlator (^ can be rewritten as 



{pip)\^{zh^[z, ~z]^{~z)\0) 



(33) 



-2impfpPf,{e* ■ z) J dye 





i{y-y)p-^ 



■'Hy) 



where h{y) = J du {tpi{u) ~ ipsiu)) , (34) 



where the ellipsis denotes terms which are responsible for 
corrections of m^-order and, therefore, they are part of 
twist 4 corrections. Q, which wc omit in this paper. 
Up to twist 3, we have 



_d_ 

dZn 



{pip)\^iz)jp[z, -z]iji-z)\0) 



= -2impfpppel j dye'^y-y^P-'hiy) . (35) 



Finally, the pure quark part of r.h.s of Eq.(pl|) together 
with the parametrization ( pi] ) and the translation invari- 
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ance yields, 



V^(-z)lO) 



-2impfpPf,e*J^ J dye 





(y) • (36) 



Thus, using (jsj,^,^), the matrix element of Eq. 
leads to 

y 1 

My) = -^f (y) - Ca^ / dx, I dx, (37) 



X2 - Xi 



y 



which reproduces the results (|27|, |25|, |2^) obtained pre- 
viously. This can be seen most easily by inserting in 
Eq.@ the expressions (H) for and ^) for 

ipf . The B terms then cancel out and what remains 
to be shown is that the r.h.s of is reproduced. The 
tpg*^" is readily obtained and the integral of ipi and of 
(fi^^ in the r.h.s of (Q) is correctly reproduced after 
interchange of integrations. This operator approach can 
be applied similarly to the axial part. 

We now relate our results (|2^, ^9|), with those obtained 
in studies of DAs based on the conformal expansion on 
the (coordinate space) light cone 0. Comparison be- 
tween (||, |l^) and similar expressions of Ref. in the 
light-cone n ■ A — gauge leads to the identification 



Biyi, y2) 
D{yi, y2) 



1 V{yi, 1 - 2/2, y2 " yi) 
mp y2 ~ yi 

1 ^(yi, 1 - y2, y2 - yi) 
rup y2 - yi 



(38) 



(39) 



Using these expressions for V and A one can now re- 
cover our results ( 27 , 2^ after straighforward but tedious 
manipulations, provided the important identification is 
done, for the vector projection: 



vi(y) = 0ii(y), f3{y) ^ g^'"Hy) 

and for the axial projection 

iag('^)(y) 



^A{y) 



dy 



(40) 



(41) 



In summary, we describe consistently exclusive processes 
at higher twist. The specific 7* — > p impact factor will 



indeed be shown to be gauge invariant and free from end- 
point singularities in a separate publication [||. The n- 
independency condition generalized up to the dynamical 
twist 3 plays a crucial role for the consistency of this 
approach with the studies of DAs performed in . 

The present framework opens the way to a system- 
atic and consistent treatment of hard exclusive processes. 
This docs not preclude the solution of the well known 
end-point singularity problem [l^ which may still re- 
quire a separate treatment. Let us finally emphasize the 
simplicity of the present method to perform practical cal- 
culations 1^. 
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